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The reason why some data on the production of if)- and T-states afeagree with 
QCD predictions, occasionally by well over one order of magnitude,!!] is that the 
traditional method for performing the perturbative calculation of the cross section 
is simply wrong. The key mistake is to require that the heavy quark pair forms a 
color singlet at short distances, given that there is an infinite time for soft gluons 
to readjust the color of the cc pair before it appears as an asymptotic tp or, alterna- 
tively, DD state. We suspect that the same mistake is made in the description of 
rapidity gaps, i.e. the production of a color-neutral quark-antiquark pair, in terms 
of the exchange of a color neutral gluon pair. The tp is after all a color neutral 
cc pair and we will show that it is produced by the same dynamics as DD pairs; 
its color happens to be bleached by soft final-state interactions. This approach to 
color is suggestive of the unorthodox prescription for the production ofrjapidity 
gaps in deep inelastic scattering, proposed by Buchmiiller and Hebecker.p When 
applied to the formation of gaps between a pair of high transverse momentum jets 
in hadron collisions, the soft color approach suggests a formation rate of gaps in 
gluon-gluon subprocesses which is similar or smaller than in quark-quark induced 
events. Formation of gaps should increase when increasing transverse momentum 
or lowering energy, in contrast with 2-gluon exchange Pomeron models. 



1 Introduction 

The conventional treatment of color, i.e., the color singlet model, has run into 
serious problems describing the data on the production of charmonium and 
upsilon statesu Specific proposals to solve the charmonium problem agree 
on the basic solution: its production is a two-step process where a heavy 
quark pair is produced first. At this stage perturbative diagrams are included 
whether the cc pair is color singlet or not. This is a departure of the textbook 
approach where only diagrams where the charm pair is in a_color singlet are 
selected. In the Bodwin-Braaten-Lepage (BBL) formalism^ the subsequent 
evolution of the pair into a colorless bound state is described by an expansion 
in powers of the relative velocity of the heavy quarks in the onium system. An 
alternative approach, color evaporation or the soft color method, represents an 
even more radical departure from the way color singlet states are conventionally 
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treated in perturbation theory. Color is, in fact, "ignored". Rather than 
explicitly imposing that the system is in a color singlet state in the short- 
distance perturbative diagrams, the appearance of color singlet asymptotic 
states depends solely on the outcome of large-distance fluctuations of quarks 
and gluons. In other words, color is a nonperturbative phenomenon. 

In Fig. |] we show typical diagrams for the production of V'-particles rep- 
resenting the competing treatments of the color quantum number. In the 
diagram of Fig. [l]a, the color singlet approach, the ip is produced in gluon- 
gluon interactions in association with a final state gluon which is required by 
color conservation. This diagram is related by crossing to the hadronic decay 
ip — > 3 gluons. In the color evaporation approach, the color singlet property 
of the ip is ignored at the perturbative stage of the calculation. The ip can, 
for instance, be produced to leading order by gg-annihilation into cc, which is 
the color-equivalent of the Drell-Yan process. This diagram is calculated per- 
turbatively; its dynamics are dictated by short-distance interactions of range 
Aa; ~ m^ 1 . It does indeed not seem logical to enforce the color singlet prop- 
erty of the ip at short distances, given that there is an infinite time for soft 
gluons to readjust the color of the cc pair before it appears as an asymptotic ip 
or, alternatively, DD state. Alternatively, it is indeed hard to imagine that a 
color singlet state formed at a range m^ 1 , automatically survives to form a ip. 
This formalism represents the original and, as we will sium. .correct method 
by which perturbative QCD calculations were performed .peffl'Q 
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(a) (b) 

Figure 1: Typical diagrams for (a) color singlet ip production and (b) color evaporation ip 
production. 



We will first discuss the resolution of the charmonium problem, emphasiz- 
ing the color evaporation approach. The solution suggests a radical departure 
from the way color is treated in perturbative QCD calculations. We will subse- 
quently speculate on the implications for the dynamics underlying the produc- 
tion of rapidity gaps which refer to regions in phase space where no hadrons 
appears as a result of the production of a color neutral quark-antiquark pair. 
The connection to charmonium physics is obvious: the ip is a color- neutral 
cc pair! 
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2 Onium Calculations with Soft Color 



The "color evaporation" or "soft color" treatment of the color quantum number 
lead to a similar description of bound and open charm production: 

1 f 2mD , da c - c 

CTonium = 7T / dm — , (1) 

9 J 2mc dm 



and 



dOrr f . da,. 



. dm — h / dm — — (2) 

9J2m c dm J 2mo dm 

9 J 2 m c dm 

where the cross section for producing heavy quarks, a cS , is computed pertur- 
batively. Diagrams are included order-by-order, irrespective of the color of the 
cc pair. The coefficients | and | represent the statistical probabilities that 
the 3x3 charm pair is asymptotically in a singlet or octet state. In order 
to achieve the phenomenological success described here it is essential to sys- 
tematically include next-to-leading order terms. Neglecting O(af) terms is 
equivalent to describing photon interactions with matter neglecting the Bcthc- 
Heitler process versus Compton scattering because it is a higher order process. 
The former actually dominates at high energy for reasons that are similar to 
those requiring the inclusion of higher order heavy quark processes. 

In principal the calculation only predicts the sum of the cross sections of 
all onium states given by Eq. (|l|). This sum rule is, unfortunately, difficult 
to test experimentally, since it requires measuring cross sections for all of the 
bound states at a given energy. This does not mean that the calculation has 
no predictive power. The above equations make the bold prediction that all 
onium states ip, tp\ \ an d r] c states share the same production dynamics which 
they also share with open charm in the limit m c ~ mu; see Eq. (3). The 
CDF collaboration has accumulated large samples of data on the production 
of prompt tp, Xcj, and tp' Since all charmonium states share the same produc- 
tion dynamics in the color evaporation scheme, their px distributions should 
be the same, up to a multiplicative constant. This prediction is borne out by 
the CDF data, as we can see in Fig. @. We will return to a detailed calculation 
of the distribution shown further on. The formalism predicts furthermore that, 
up to color and normalization factors, the energy, xp and px dependence of 
production cross sections for onium states and open charm pairs is the same. 
Support for the prediction of Eqs. ([!]) and (||) that the production of hid- 
den and open charm have similar dynamics is shown in Fig. |^, which displays 
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Figure 2: Data from the CDF Collaboration, H shown with arbitrary normalization. The 
curves are the predictions of the color evaporation model at tree level, also shown with 
arbitrary normalization. The normalization is correctly predicted within a K-factor of 2.2. 



charm photoproduction data for both open charm and bound state production 
with common normalization in order to show their identical energy behavior. 
A similar figure for hadroproduction can be found in Ref. 11. By the same 
argument the formalism also predicts that the normalized xp distribution for 
J/ip and DD paiKS-shpuld be the same at a given center-of-mass energy. This 
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Figure 3: Photo-production datapB as a function of the photon energy in the hadron rest 
frame, W 1 . The normalization has been adjusted to show the similar shapes of the data. 



One of the most striking features of color evaporation is that the produc- 
tion of charmonium is dominated by the conversion of a colored gluon into 
a "0, as in Fig. [l]b. In the conventional treatment, where color singlet states 
are formed at the perturbative level, 3 gluons (or 2 gluons and a photon) are 
required to produce a ip. Contrary to the usual folklore, "0's are, except at the 
higher energies, not produced by gluons. As a consequence color evaporation 
predicts an enhanced ip cross section for antiproton beams, while the color 
singlet model predicts roughly equal cross sections for proton and antiproton 
beams. The prediction of an enhanced p yield is obviously correct: antiproton 
production of ip's exceeds that by protons by a %jtftr 5 close to threshold; see 
Fig. [| This fact has been known for some timeBErtl We should note that for 
sufficiently high energies, gluon initial states will eventually dominate because 
they represent the bulk of soft partons. 

3 Quantitative Tests of Soft Color 

The color evaporation scheme assumes a factorization of the production of the 
cc pair, which is perturbative and process dependent, and the materialization 
of this pair into a charmonium state by a mechanism that is nonperturbative 
and process independent. This assumption is reasonable since the character- 
istic time scales of the two processes are very different: the time scale for the 
production of the pair is the inverse of the heavy quark mass, while the for- 
mation of the bound state is longer than the time scale 1/Aqcd- Therefore, 
explicit comparison with the tp data requires knowledge of the fraction of 
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Figure 4: Normalized xp distribution for the production of of DD (histogram) at y/s = 27.4 
GeV and J/ip (squares) in proton- nucleon collisions at y/s = 23.7 GeV. The curves are the 
prediction of our model. 



produced onium states that materialize as ?/>'s, i.e., 



where p^ is assumed to be a constant, independent Jif the process. This as- 
sumption is in agreement with the low energy dataJi3 : t£l The constant not only 
accounts for the direct production of ip to the onium cross section, but also 
includes its production via ip' and x production and decay. 

Quantitative tests of color evaporation are made possible by the fact that 
all ^-production data, i.e. photo-, hadroproduction, Z-decay, etc., are de- 
scribed in terms of a single parameter. Once has been empirically de- 
termined for one initial state, the cross section is predicted without free pa- 
rameters for the other. We will illustrate the power of the color evaporation 
scheme by showing how it quantitatively accommodates all measurements, in- 
cluding the high energy Tevatron and HERA data, which have represented a 
considerable challenge for the color singlet model. Its parameter-free predic- 
tion for the Z-boson decay rate into ip's is an order of magnitude larger than 
the color singlet model and consistent with data. 

In Fig. U we compare the photoproduction data with theory, using the NLO 
perturbative QCD calculation of charm pair production from Ref. 17. From 
the relative magnitude of the ip and open charm cross sections we determine 
the fragmentation factor p^ to be 0.50 using GRV HO, or 0.43 using MRS A 
structure functions. Note that the factor p^ possesses a theoretical uncertainty 
due to the choice of scales and parton distribution functions. We conclude the 
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Figure 5: Ratio of the cross sections for the production of J/ip by proton and antiproton 
beams in the color evaporation model (solid line) and the color singlet model (dashed line) 
as a function of the center-of-mass energy. Data taken from Ref. 14. 



photoproduction of J/ip and DD is well described by the color evaporation 
model. This reaction has now been used to fix the only free parameter, ps 
0.5. 

At this point the predictions of the color evaporation model for hadropro- 
duction of ip are completely determined, up to O(ctg) QCD corrections. In 
Fig. |?] we compare the color evaporation model predictions with the data, and 
conclude that the this color scheme describes the hadroproduction very accu- 
rately. In order to also obtain a theoretical prediction for the pr-distribution 
already shown in Fig. ||, we have computed the processes g + g — > [cc] + xu 
q + q — » [cc] + g, and g + q — ► [cc] + q at tree level using MADGRAPH.I13 
We imposed that the cc pair satisfy the invariant mass constraint of Eq. ([!]). 
Our results are shown in Fig. p[ Higher order corrections such as soft-gluon 
resummation are expected to tiLL our lowest order prediction, bringing it to a 
closer agreement with the data.t3 

In the color-evaporation scheme the width for inclusive Z decay into prompt 
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Figure 6: Photoproduction dataotiJ and the predictions of the color evaporation model at 
next-to-leading order as a function of the photon energy in the hadron rest frame, Wy. The 
normalizations in this figure are absolute. 



charmonium is: 

1 f 2m ° dt ce 
T(Z~^ prompt charmonium) — — dm —r— , (5) 

9 J2m c dm 

where T is the partonic width for producing a cc pair. The procedure should by 
now be familiar: in order to obtain the partial width into a specific charmonium 
state we multiply the above expression by the appropriate fragmentation frac- 
tion p into ip, which was determined from charmonium photoproduction data. 
Notice that the predictions for the Z decay into charmonium are parameter- 
free. 

We have again evaluated all the tree-level partonic amplitudes using the 
package MADGRAPH.E3 Although formally of higher order in a s , the domi- 
nant process for the inclusive decay of the Z into charmonium is Z — > ccqq, 
where q = u, d, s, c, and b. (The leading-order process in a s is Z — > ccg, which 
leads to the production of a charmonium state and a hard jet is suppressed by 
a virtual quark propagator of order m c /mz). The branching fraction of Z into 
prompt ip is (1.7-1.8) x 10~ 4 . This is to be contrasted with the color-singlet 
model which predicts a branching fraction for direct ij) in Z decay of the order 
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Figure 7: Hadroproduction dataEj'til and the predictions of the color evaporation model 
at next-to-leading order as a function of the center-of-mass energy, E cm - The curve for 
bound state production is an absolutely normalized, parameter-free prediction of the color 
evaporation model. 
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3 x 10~ 5 .E3 The color-evaporation model leads to a branching fraction larger 
by almost an order of magnitude consistent with the result reported by the 
OPAL collaboration of 

B(Z -> prompt ^ + X) = (1.9 ± 0.7 ± 0.5 ± 0.5) x 10~ 4 . 

We hope that we have illustrated by now that the soft color approach gives 
a complete picture of charmonium production in hadron-hadron, 7-hadron, and 
Z decays. The phenomenological success of the soft color scheme is impressive 
and extends to applications to other charmonium and upsilon statesJij'tj 

4 Intermezzo: Soft Color and BBL 

Other approaches, very similar in spirit, can be found in Refs. 3, 15 and 23. 
The color evaporation approach differs from Ref. 3, the formalism of Bodwin, 
Braaten and Lepage, in the way that the cc pair exchanges color with the 
underlying event. In the BBL formalism, multiple gluon interactions with 
the cc pair are suppressed by powers of v, the relative velocity of the heavy 
quarks within the tp. The color evaporation model assumes that these low- 
energy interactions can take place through multiple, soft-gluon interactions. 
While the formalism allows straightforward application to heavy quark decays, 
it is not always clear how to compute production cipss sections in the BBL 
formalism (e.g. photoproduction of ip near z = l)Ji3 The color evaporation 
scheme, though partly nonperturbative, is phcnomcnologically well-defined, 
has less parameters (1 versus 3 for describing i/'-production). Also, next-to- 
leading order corrections are included in a straightforward way, a necessary 
condition for obtaining quantitative predictions. 

The ^>'s produced through the color-evaporation mechanism are expected 
to be unpolarized since the polarization information is lost because of the 
multiple soft gluon exchangesO On the other hand, the (nonlppLarizatiorj of 
ip is hard to explain in the framework of the color-octet model.Ejo Therefore, 
the measurement of the polarization of the produced charmonium may very 
well be a tool to discriminate between these competing descriptions. 

5 Implications for the Physics of Rapidity Gaps 

The important lesson about color resides however in the similarity, not the 
differences of these approaches: perturbative color octet states fully contribute 
to the asymptotic production of color singlet states such as tp's. We sus- 
pect that this is also true for the production of a rapidity gap which is, e.g. 
when produced in electroproduction, nothing but the creation of a color singlet 
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Figure 8: Mechanism for the production of rapidity gaps in deep inelastic scattering. 



quark-antiquark pair; see Fig. ^. The diagram shown represents the produc- 
tion of final state hadrons which are ordered in rapidity. From top to bottom 
we find the fragments of the intermediate partonic quark-antiquark state and 
those of the target. Buchmullcr and Hebecker proposed that the origin of a 
rapidity gap corresponds to the absence of color between photon and proton, 
i.e. the 3 x 3 (= 1 + 8) intermediate quark-antiquark state is in a color singlet 
state. Because color is the source of hadrons, only the color octet states yield 
hadronic asymptotic states. This leads to the approximate expectation that 

F^ P) = j^F, (6) 

Although this result is subject to corrections, it embodies the essential physics: 
events with and without gaps are described by the same short-distance dynam- 
ics. Essentially non-perturbative final-state interactions dictate the appearance 
of gaps whose frequency is determined by simple counting. The treatment of 
color is the same as in the case of heavy quark production: the same pertur- 
bative mechanisms, i.e. gluon exchange, dictates the dynamics of color-singlet 
gap (ip) and regular deep inelastic (open charm) events. 

Our understanding of the (soft) nature of color challanges the orthodox 
description of rapidity gaps in terms of the so-called hard Pomeron description 
sketched in Fig. [| The i-channel exchange of a pair of gluons in a color 
singlet state is the origin of the gap. The color string which connects photon 
and proton in diagrams such as the one in Fig. pi is absent and no hadrons 
are produced in the rapidity region separating them. The same mechanism 
predicts rapidity gaps between a pair of jets produced in hadronic collisions; 



see Fig. 10. These have been observed and occur with a frequency of order of 
one percent. Ej The arguments developed in this work invalidate this approach: 
it is as meaningless to enforce the color singlet nature of the gluon pair as 
it is to require that the cc pair producing ip is colorless at the perturbative 
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Figure 9: Pomeron mechanism for the formation of rapidity gaps. 

level. Following our color scheme the gaps are accommodated as a mere final 
state color bleaching phenomenon a la Buchmiillcr and Hcbecker. This can be 
visualized using the diagram shown in Fig. |ll[ At short distances it represents 
a conventional perturbative diagram for the production of a pair of jets. Also 
shown is the string picture for the formation of the final state hadrons. Color 
in the final state is bleached by strings connecting the 3 jet at the top with 
the 3 spectator di-quark at the bottom and vice-versa. The probability to 
form a gap can be counted a la Buchmiillcr and Hebecker to be 1/(1 + 8) 2 
because it requires the formation of singlets in 2 strings. This is consistent with 
observation and predicts that, as was the case for electroproduction, the same 
short-distance dynamics governs events with and without rapidity gaps. The 
datanJ is consistent with the prediction of this simple picture which basically 
predicts that the gap fraction between pp jets is the square of that between 
virtual photon and proton in deep inelastic scattering. 

, jet 




Figure 10: Pomeron mechanism for the formation of rapidity gaps in hadron collisions. 

One should realize that this string picture is not necessarily the correct 
one. It is more likely that the color is bleached between the top and bottom 3 
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Figure 11: Color bleaching picture for the formation of rapidity gaps in hadron collisions. 



and 3 which are widely separated in rapidity space. 

This discussion ignores that gluon-gluon as well as quark-quark subpro- 
cesses contribute to jet production in hadron collisions. In the color flow dia- 
gram corresponding to Fig. [ll] top and bottom protons each split into a color 
octet gluon and color octet 3-quark remnant. There are now (8 x 8) 2 color 
final states. We anticipate a reduced probability to form a color singlet. The 
reduction may not be very significant because the 10 + 10 and 27 color final 
states may be suppressed. One argument for this is that these representations 
consist of exotic multi-quark states which do not materialize into final state 
mesons. High color charges may also be suppressed for dynamical reasons. 
Despite the fact that we can at best guess the non-perturbative dynamics, it 
is clear that the soft color formalism predicts a gap rate which is similar of 
smaller in gluon-gluon interactions. This is in contrast with the diagram of 
Fig. ^ which predicts a gap-rate enhanced by a factor (|) in gluon-gluon sub- 
processes.ES The contrasting predictions can be easily tested by enhancing the 
relative importance of quark-quark subprocesses. i.e. by increasing the pr of 
the jets at fixed energy of by decreasing the collision energy of the hadrons 
at fixed px- In either case we anticipate in the soft color scheme an increased 
rate for the production of gaps, a prediction opposite from that expected in 
the 2-gluon exchange model. o 

Do Figs. [§l-|lT| suggest that we have formulated alternative s- and i-channel 
pictures to view the same physics? Although they seem at first radically dif- 
ferent, this may not be the case. Computation of the exchange of a pair 
of colorless gluons in the i-channel is not straightforward and embodies all 
the unsolved mysteries of constructing the "Pomeron" in QCD. In a class of 
models where tjac-Bomeron is constructed out of gluons with a dynamically 
generated masspJ'H the diagram of Fig. ^ is, not surprisingly, dominated by 
the configuration where one gluon is hard and the other soft. The diagram 
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is identical to the standard perturbative diagram except for the presence of a 
soft, long-wavelength gluon whose only role is to bleach color. Its dynamical 
role is minimal, events with gaps are not really different from events without 
them. Soft gluons readjust the color at large distances and long times. Their 
description is outside the realm of perturbative QCD. In this class of models 
the hard Pomeron is expected to be no more than an order ct2 correction, a 
view which can be defended on more solid theoretical ground.E3 

Some have challenged the theoretical soundness of this line of thinking.OO 
Also note that our discussion is at best indirectly relevant to completely non- 
perturbative phenomena like elastic scattering. There is no short distance limit 
defined by a large scale. The Pomeron exists. 
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